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Abstract: 

We apply the optimization procedure based on the Principle of Minimal Sensitivity to 
the third-order calculation of R-t-. Since the effective couplant remains finite, freezing to a 
value Os/vr = 0.26 at low energies, we can actually evaluate the defining integral of Rr and 
compare the optimized perturbation theory result to that of the optimized result obtained 
after the integral has been evaluated using contour techniques. The good agreement shows 
that the optimization procedure is consistent and suggests that the infrared fixed point is 
meaningful. 



1 Introduction 

The inclusive semihadronic decay rate of the r lepton, expressed as the ratio 

r(r" i^r + hadrons) 

is a fundamental test of QCD. The parton model gives a rough estimate by approximating 
hadronic decay as the rate into quark-antiquark pairs: Rr ^ Nc = 3. Corrections to this 
approximation include both perturbative and non-perturbative QCD, as well as electroweak 
corrections |2[. Here we will concentrate on the pertubative QCD corrections and compare 
two different ways of applying the optimized perturbation theory to Rr- 

Refs. m, |2| discuss in detail the theoretical calculation of Rr, which has the following 
form: 

1^ = 1 ^sD(s)^l <i»^(l-jj,) (l + ^)fiW. (2) 

with 

R{S) = 3(|Kd|' + |Ks|')[l +^r(s)]. (3) 

Using the couplant a = ^ the perturbative corrections are expressed as a series, truncated 
to third-order: 

= a(l + ria + r2a^). (4) 

From the recently re-done third-order calculation of Gorishny, et al. |^, Q : 

ri(MS,;U = M^) = 1.986 - 0.115n/, (5) 

r2(MS,;U = Mr) = -6.637 - 1.200n/ - 0.005n}, (6) 

in the MS scheme with the renormalization scale /x taken to be Mr- For three flavors, 
ri = 1.64 and r2 = —10.28. [Note this differs greatly from the earlier result r2 = 93.98 |l[| 
based on the erroneous calculation of |^].] 
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Historically the integral @ has not been evaluated directly because it was expected 
that the QCD couplant would become large at small s, making a perturbative prediction 
impossible. Instead, Rr was re-expressed as a contour integral in the complex s plane with 
the contour running clockwise around the circle of |s| = M^, thus avoiding the small-s 
region. To evaluate the integral, the coupling constant is expanded in powers of asiM-j-) 
The final result (for three flavors) is then, in the MS(// = M,-) scheme, 

Rr = 3{\Vudf + \Vusf)[l + a + 5.20a2 + 26.37a=^]. (7) 
2 Optimized Perturbation Theory 

Based on the principle of minimal sensitivity, optimized perturbation theory (OPT) Q finds 
the renormalization scheme (RS) in which the result is least sensitive to changes in the RS 
parameters. For a detailed discussion of the OPT method see The application of Ref. 
10] to the Re+e- ratio in third order is easily adapted to the present case. In effect, the 
difference between QCD corrections to the ratio of the e~^e~ hadronic cross section and 
the r hadronic decay involves taking (J^Qf)'^ — (i^ f2{MS)) and replacing 3X]Q/ by 
3EV7/«3. 

As previously explained in |p the question of an infrared fixed point can be addressed 
by solving the optimization equations and requires a third-order calculation to determine 
the RS invariant p2-, where 

P2 = r2+C2- (ri + -cf, (8) 
and c and C2 are /3-function coefficients; 

IJi— = (3{a) = -ha^{l + ca + 020^ + ...). (9) 
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If p2 is negative a positive fixed point, a*, exists, and the more negative p2 is, the smaher 
that a* will be. Just as in the R^+e- case, here p2 is negative and an infrared fixed point 
exists. Our analysis shows that the numerical solutions to the optimization equations tend 
towards that fixed-point solution. Because of this "freezing" of the coupling constant to a 
small finite value, the perturbative expansion of 7^t-(s) does not diverge and the integral (|2|) 
can be evaluated numerically. Thus there are two ways of applying optimized perturbation 
theory to the problem of evaluating Rr '- (i) optimization of the integrated expression for i?^, 
i.e. Eq. (|^, and (ii) optimization of the perturbative expansion inside the integral down 
to low energies and then integrating the result. The second method gives a perturbative 
prediction for the differential decay rate. 

Method (i) involves applying optimized perturbation theory to the MS result for i?,-, (j^, 
and was previously done by Chyla et ai, |Q. Here we take Vud = 0.9747 and Vus = 0.218 
from the particle data book [^]. For three flavors and s = [using AgL = 280 MeV 
and Mr = 1.777 GeV we get the following optimum values: a^^^Mr) = 0.164, 
fi = -0.49, f2 = 2.98 and Rr = 3.48. 

Now we can compare this result to method (ii) where we have to solve the optimization 
equations for TZris) in the range > s > and evaluate the integrand of Eq. (^) explicitly. 
Note that the perturbative expression to be optimized in this case is (^) which has different 
coefficients than (0). At s = the relevant number of flavors is three (and again we will 
take A^ = 280 MeV for comparison to method (i)) but at lower energies there are only two 
active flavors, so we must 'match' to a 2-flavor theory with an appropriate A. We follow the 
method outlined in by adjusting ^jj^ so that TZr is continuous at the matching point, 
which we take to be ^/s = 2ms, (with a current quark mass of 199 MeV it 33 [^]). This 
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leads to a value of = 250 MeV. Another relevant point is that the term |VusP must 
be set to zero below the s,u kinematic threshold at ^/s = nis + m^, since for ^/s below this 
value the virtual W can decay only to an u,d pair. 

As with i?g+g_, solving the optimization equations is straightforward except for the slow 
convergence of the iteration method at low energies. Numerical solutions were obtained 
down to ^/s = 0.10 GeV, where the result joins smoothly to the analytic result for the 
^/s — > limit 0]. In Fig. 1 we show the optimized couplant a as a function of ^/s. Also 
shown is the optimized third-order correction 1Zr{s). At zero energy we have the infrared- 
fixed point solution for two flavors. In Table 1 we list the infrared values for 2 and 3 flavors. 
[In the Re+e- case, the results p2 and a* are identical for 3 flavors (since i^Qf) = for 
u^d,s quarks), but slightly different for 2 flavors {p2 = —10.91, a* = 0.263).] 

From the optimized result for Rr{s) we can obtain D{s) from Eqs. (§),(|3|). This is 
plotted in Fig. 2. Integrating this from to M2 gives Rr = 3.44, which agrees well with 
the result 3.48 obtained with method (i). 

D{s) represents a perturbative prediction for the differential decay rate of dV{T~ — > 
v-r -|- hadrons(s))/(is, normalized by r(r — > VT-e~iJe), for r decays to hadrons with an invari- 



ant mass of ^/s [13|. We do not expect this prediction to be right, of course. The corre- 
sponding experimental quantity would presumably show structure due to meson thresholds 
and resonances, particularly vr and ai. However our hypothesis is that the perturbative 
prediction is meaningful in that if both the data and the prediction are 'smoothed' in some 
suitable fashion then they will agree. This hypothesis is supported by the results of Ref.0] 
which finds excellent agreement when Poggio-Quinn- Weinberg (PQW) smearing is applied 
to both the experimental data and the theoretical prediction for -Rg+e-- We suggest that 
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PQW smearing applied to R{s), related to D{s) by (|2|), will also give good agreement in 
the r case. Unfortunately, it does not seem possible to test this prediction at present, since 
existing r-decay data give inadequate information about neutral hadrons. It is a bold predic- 
tion because it implies that the vector and axial-vector contributions, although dominated 
by different hadrons, must become the same after smearing. While the smearing is quite 
drastic (since it must smooth out the peak structure in the data) the theory /experiment 
comparison it allows is quantitative and highly non-trivial. 

3 Theoretical Uncertainties 

To answer the question of how good the agreement is between the two methods we need to 

(3) 

ask where the uncertainty lies in the calculation. Other than the choice of (which was 
the same for both methods) there is only one uncertainty in method (i): the truncation of 
the perturbative series. Following the argument of Q we estimate the error as |f2a^|, and 
thus for method (i): 

= 3.48 ± 0.04. (10) 

In method (ii) there are two sources of error: truncation of the perturbation series 
and uncertainty in the input parameters (the strange quark mass and A^^^). The series- 
truncation error was again estimated as \f2a^\ and then integrated as in Eq. (^). This had 
an uncertainty of about ±0.10. 

As explained by Marciano [0] the A^ of our effective theories (massless quarks with 
different numbers of flavors) must be matched at the thresholds to correspond to a single, 
underlying "full QCD" theory. The exact procedure for doing this, though, is unknown; 
should one make TZ continuous or continuous across the threshold? The differences 
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between the two procedures are small, only about 5 MeV, but we make the following estimate 

(2) 

to gauge the size of error: if A^g^ increases to 260 MeV then Rr increases by only 0.002, 

(2) 

which shows that the perturbative corrections are insensitive to this uncertainty in ^j^s- 
This agrees with the e~^e~ analysis, 0, which found that below the s-quark threshold Rg+g- 
was most sensitive to the fixed-point solution. 

The estimated error on the current strange quark's mass is about 20%. If the mass of 
the strange quark is decreased to 166 MeV, the 3 flavor region now runs down to 332 MeV. 
The effects of this on TZr are barely discernable and do not significantly modify the integral. 

(2) 

The moving of the threshold also affects the value of AV^. For both limits of the strange 

MS ^ 

(2) 

quark mass the A^-rrr increases to about 260 MeV. As we found above this also has little 

^ MS 

effect on the error. The net effect of changing the strange quark mass is estimated to be 
±0.006. Therefore with method (ii) we estimate: 

= 3.44 ±0.11. (11) 

As one can see the two methods agree well within the error estimates of method (ii). 

4 Comparison to Experiment 

Rr can be found by measuring the leptonic branching fractions, Be, B^^ Q 
Independently Rr can also be found by measuring the total decay rate: 



i^r '^T '^T^e Tt-_>^ 



r 

T— >e 

where T^-^i = r(T~ Ut-1~i>i) can be calculated theoretically with great accuracy because 
it is a purely electroweak process. Using a mean life time of l/F = tq = (0.3025 ± 0.0059) x 
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IQ-^^s Ref. § quotes i?f ^'^ = 3.32 ± 0.12. Taking the average of these two values gives: 

RfP = 3.61 ± 0.05. (13) 

To compare to experiment we need to estimate the size of error due to the choice of 
AVr;-.. Increasing the value by 50 MeV to 330 MeV leads to a 2% increase in method (i) and 
1% increase in method (ii). In Table 2 we list the actual values for the case of increasing 
and decreasing A^g^ by 50 MeV. Note that method (i) is more sensitive to these changes 
than method (ii). 

It is generally accepted that non-pertubative corrections are small, decreasing R by only 
0.5% RP'^^^ [|l5|. Electroweak corrections are estimated to increase R by 2%. Therefore we 
get the following estimates for the two methods: 

R^;^ = 3.48 - 0.02 + 0.07 = 3.53 ± 0.12, (14) 

jliu) ^ 3 44 _ Q_Q2 ^ 0.07 = 3.48 ± 0.15. (15) 

Both methods give very good results, which are near the experimental average and well 
within the errors. 

5 Conclusion 

Using optimized perturbation theory we have shown that the infrared fixed point and the 
'freezing' of the strong coupling constant at low energies can be used to evaluate the Rr 
integral directly. Comparison to the usual OPT method of optimizing the MS result is 
excellent showing the consistency of OPT. Futhermore both methods give good agreement 
to the experimental value. 
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nj f>2 a* UriO) 



2 -10.83 0.264 0.332 

3 -12.21 0.244 0.303 



Table 1: Infrared fixed point values for 2 and 3 flavors. 



(MeV) 


7^i') 


7^i") 


330 


3.56 ± 0.06 


3.48 ± 0.13 


280 


3.48 ± 0.04 


3.44 ± 0.10 


230 


3.41 ± 0.02 


3.40 ± 0.08 



Table 2: Comparison of methods (i) and (ii) for different ^j^- Errors are estimated 
by |r2a^|, integrated for method (ii). 
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Fig. 1. The optimized third-order results for a = Os/tt and TZj- {s). The vertical line 
indicates where the Nf = 2 and Nf = 3 effective theories are matched. Error bars for 

1Z^\s) are ±|f2a^|. 



Fig. 2. Integrand of Rr; i.e. the differential decay rate dr(T~ 1/-^ + hadrons(s))/ds for 
decays into hadrons with invariant mass-squared s, divided by r(r~ — > i'r^~i'e)- The 'jump' 
at -s/s 200 MeV occurs because the term \Vus\^ = below the s,u kinematic threshold. 
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